Coulomb interaction, ripples, and the minimal conductivity of graphene 



00 
O 

o 

(N 

oo 
(N 



CD 



o 

> 

o 
o 



X 




Igor F. Herbut^, Vladimir Juricic^, and Oskar Vafek^ 
^ Department of Physics, Simon Eraser University, Burnaby, British Columbia, Canada V5A 1S6 
^ National High Magnetic Field Laboratory and Department of Physics, 
Florida State University, Tallahassee, Florida 32306, USA 

We argue that the unscreened Coulomb interaction in graphene provides a positive, universal, and 
logarithmic correction to scaling of zero-temperature conductivity with frequency. The combined 
effect of the disorder due to wrinkling of the graphene sheet and the long range electron-electron 
interactions is a finite positive contribution to the dc conductivity. This contribution is disorder 
strength dependent and thus non-universal. The low-energy behavior of such a system is governed 
by the line of fixed points at which both the interaction and disorder are finite, and the density of 
states is exactly linear. An estimate of the typical random vector potential representing ripples in 
graphene brings the theoretical value of the minimal conductivity into the vicinity of 4e^ /h. 
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Graphene, an atom-thick layer of graphite, has re- 
cently defined a new frontier of condensed matter physics. 
Its essential electronic property, inherent to the honey- 
comb lattice formed by the carbon atoms, is that the 
low-energy quasiparticle excitations can be thought of as 
being massless Dirac fermions, which propagate with a 
Fermi velocity of around three thousandths of the ve- 
locity of light. This pseudo-relativistic nature of the 
quasiparticle excitations makes the electronic proper- 
ties of graphene fundamentally new in many respects 
[ij. When the chemical potential is tuned to the Dirac 
point graphene provides a rare example of a critical 
two-dimensional fermionic system In contrast to its 
textbook bosonic equivalent all sufficiently weak in- 
teractions between electrons, including the long^anged 
Coulomb, are then irrelevant perturbations IHl. The 
effects of electron interactions thus become progressively 
less important as the system is probed at lower frequen- 
cies and temperatures. One important consequence of 
this "infrared freedom" is that the zero-temperature dc 
conductivity of clean graphene is finite and universal, and 
simply determined by its gaussian value of = (7r/2)e^//i 
[il . Including scattering of impurities in a self-consistent 
Born approach yields another, similar in magnitude and 
still universal, value of (4/7r)e^//i [2], [J. Localiza- 
tion corrections are also expected to set in at very low 
temperatures [8] and thus further diminish the conduc- 
tivity. The actual measurements of graphene's conduc- 
tivity, however, are in significant discord with these re- 
sults: experimentally, a ~ Ae'^/h, and thus significantly 
larger than all the theoretical values The origin of 
this discrepancy is unclear at the moment, with several 
recent works focusing on the role played by the extrinsic 
charged impurities [9|. 

In this Letter we show that the long-range Coulomb 
interaction between electrons in graphene provides the 
leading correction to the gaussian value of conductivity, 
which is positive and by itself only logarithmically slowly 
vanishing when frequency approaches zero. This suggests 



FIG. 1: The flow of the disorder strength A characterizing 
ripples of the graphene sheet and the Coulomb coupling g at 
weak couplings. Two effects are balanced at the attractive line 
of fixed points that emanates from the gaussian fixed point 
at the origin. At the line the system exhibits linear density 
of states and a finite non-universal dc conductivity. 



that the origin of the observed unusually large conduc- 
tivity at the Dirac point may be intrinsic, and originate 
from the Coulomb correction which is then effectively cut 
off by a finite temperature/disorder/size effects. 

We consider a specific mechanism of such cutoff which 
invokes a non-trivial int erp lay between long-range in- 
teractions and disorder j^Oj, the latter being presently 
provided by the apparently unavoidable wrinkling of the 
graphene sheet. The combination of the Coulomb inter- 
action and the random vector potential that may be used 
to represent such ripples in graphene leads to a line of sta- 
ble fixed points (see Fig. 1). The finite zero-temperature 
dc conductivity, obtained here from the Kubo formula, 
varies continuously along the line and, most importantly, 
increases with the increasing disorder strength. We pro- 
vide the symmetry arguments for the existence of the 
line of stable fixed points, and extract the density of 
states at low energies. A crude estimate of the typical 
parameters in graphene gives a sizable correction to the 
gaussian value already to the lowest order in our calcu- 
lation and significantly narrows the gap that presently 



exists between the theory and the experiment. Further 
impUcations of our theory and the connections with the 
critical bosonic theories and related theoretical results in 
literature are discussed. 

The low-energy excitations in the vicinity of the Dirac 
points at aK, ct = ±1 may be represented by two four- 
component Dirac spinors = (^'j^^(x, r), *I/^^(x, t)). 
^ctT is a two-component Grassman field representing the 
components of the electron with the third component of 
its (real) spin up on the two sublattices of the honey- 
comb lattice and with wavevectors near aK The 



imaginary-time Lagrangian density for the interacting 
system of quasiparticles in presence of a random vector 
potential representing ripples of the graphene sheet is 
then L = Lq + Lc + Ld, where 



and 



(1) 



(2) 



(3) 



Here the index a = 1,2..., TV labels replicas introduced 
to average over disorder, = 0,1,2, and n = 1,2, and 
{Tmi 71^} = 2(5^j,. The summation over repeated indices is 
assumed and the limit iV ^ is to be taken at the end. 
The integration over ao (x, t) reproduces the standard ^ 
(7/|x| electron-electron interaction, whereas integrating 
out the static vector potential An (x) yields an alternative 
form oi Ld [6]: 

^D^jJ dr'(a*^7n*^)(x,T)(a'*f,7n*f,)(x,T'), 

(4) 

which will be also used. For convenience, we have set h — 
e/c = vp — 1, where vp « lO^m/s is the Fermi velocity. 
In this convention there are two dimensionless coupling 
constants in the theory: the Coulomb interaction g = 
27re'^ /ehvp, and the strength of disorder A. 

In analogy with the two-dimensional bosonic critical 
systems Q , at T = the conductivity at frequencies well 
below the microscopic energy scale ft/b in the units of 

/h can be written in the scaling form as 



a{Lo)^F{bcu,g{b),A{b)), 



(5) 



where F(x, y, z) is a universal scaling function. The func- 
tional dependence oi g{b) and A (6) ensures the ultimate 
independence of (j{ui) on the arbitrary factor b. We have 
set all the irrelevant couplings to zero. Let us choose 
then bu) = K <^Vl, with A as an arbitrary scale. Then 



FIG. 2: The one- loop corrections to the disorder vertex in 
Eq. (4). The dashed and wiggly lines stand for the scalar 
and the vector field propagators, respectively. The last three 
diagrams sum up to zero. 



If the couplings g{A/uj) and A(A/a;) flow to small values 
as a; ^ 0, one can expand the function F{A, g, A) as 

a{uj) = F(A, 0, 0) + ug (A/lo) + wA [A/uj) , (7) 

to the first order in the two coupling constants, u and 
w are constants. Kubo formula yields the value of the 
gaussian term F(A, 0, 0) = 7r/2 in the continuum limit @ 
(see also Eq. (16) below). In what follows we compute 
the remaining two terms in the last expression. 

Let us first obtain the cutoff-dependent couphngs g(b) 
and A(&). Non-analitycity in momentum of the second 
term in Lc implies that the flow of g with the change of 
cutoff can be written exactly [lOj as. 



dg 



d\n{b) 



= (^-1)5, 



(8) 



with the dynamical exponent z fixed by the requirement 
b'^~^ = Zi^/Zk- Z^ and are the wave-function and the 
velocity renormalizations, respectively. To the first order 
in both couplings one finds 
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(9) 



in agrement with previous calculations f^, [l^ . 
[l^. At g = 0, the disorder strength A is exactly 
marginal coupling Q . We have also confirmed this by an 
explicit calculation to order A^ (Fig. 2). When g ^ 0, 
A (6) = AZ/^/Zl. Using the form in Eq. (4) and com- 
puting the first diagram in Fig. 2 then gives 



ZA = l + i^ln(6). 
47r 



(10) 



Thus, as the ultraviolet cutoff in the theory is changed 
from r2 to 11/6, the Coulomb and the disorder couplings 
flow according to the differential equations [l2| 



J|5y^»(^^+o(.^A^.A,|, (11, 



dA 
dln(6) 



= 0. 



(12) 



a{Lo) = F{A,g{A/u,),A{A/Lo)). 



(6) 



Under renormalization the electron interaction may thus 
both decrease or increase, depending on disorder. 

Although the calculation has been performed here only 
to the leading order, we suspect that the equation (E 



3 



may in fact be exact. Without the last term in Lu the 
rest of the (interacting) Lagrangian L enjoys the symme- 
try under the time-independent gauge transformation 

A„(x)^A„(x)+a„x(x), vl;^ ^ e-x(x),j,a_ (^3) 

This imphes that the usual Ward identities hold and that 
the polarization of the vector field A„(x) is transverse. 
The coupling A, which in Ljj appears as the inverse mass 
for A„ 



i(x), should therefore not renormalize 
For completeness let us note that L is also symmetric 
under the purely time-dependent gauge transformation 



af^(x,T) 



a^{^,T)+drr{T), (14) 



which guarantees the preservation of the form of Lc and 
is ultimately responsible for Eq. (8) [3]. 

The conductivity to the first order in g and A may 
be computed next. We couple the external electromag- 
netic vector potential a minimally to Dirac fermions and 
choose a = aei, for example. Using the Kubo formula 
for the replicated theory [15| and to the first order in the 
two couplings we find 

F{-iLu,g,A) ^ Ig{lo) + glciio) + AId{lo), (15) 

where the gaussian value is 

Ig{uj) = 167r— / .„ ^„ — ^ (16) 



dujj (27r)3 q^{{qo-ujf + ce) 



= _+0(_ 

2 W 



To the first order in disorder only the self-energy diagram 
contributes and yields a negative contribution 



^ , , „ . d f dvd 'k 

/.(.)^64.-y^ 



dvd'^kd^p v{v'^ — 



f (iy2 +k2)(zy2 +p2)2 



(17) 



Finally, the Coulomb contribution is 
d f d^kd^p 



Ic{uj 



1 /t^ 

- + o{- 

6 \n 



4k - p 



4k2 uj^ + 4p' 



duj J (27r)3 |k + p|(tj2 _^4k2 
2 ' 4k • pi /25 7r\ 

48 ^ 8 



18) 



Few comments are in order at this point. The deriva- 
tives with respect to frequency in Eqs. (16)-(18), which 
are to be taken for strictly positive frequencies, serve to 
subtract the finite uj — contributions to the integrals. 
These are known to arise as the artifacts of the ultravi- 
olet cutoff, which violates gauge invariance. Second, in 
contrast to the disorder contribution, the Coulomb term 
contains both the self-energy and the vertex corrections, 
given by the first and the second term in the curly bracket 



in Eq. (18), respectively. While each of these two sepa- 
rately is logarithmically divergent in the continuum limit 

— > 00, the divergences cancel out exactly in the full ex- 
pression for the conductivity. The final result in Eq. (18) 
represents the finite remnant left after this cancelation. 
Indeed, the cancelation of logarithms in the conductiv- 
ity is to be expected: without it the field-theoretic result 
would not be cutoff-independent, and would be physically 
meaningless. The results in Eqs. (16)-(18) are universal 
numbers characteristic of the continuum limit. 

The last three expressions together with Eqs. (7), (11), 
and (12) at finite disorder then give the dc conductivity 



a(0) = 



(4 - 7r)A + 0(A2) 



h ' 



(19) 



and non-universal. The result is ^an^er than the Gaussian 
value due to the positive Coulomb contribution at the line 
of fixed points. Note that A in the last expression is the 
same as the "bare" value of the disorder strength at the 
microscopic scale. 

In the ideal sample with A = 0, on the other hand, 
solving the Eq. (11) and inserting into Eq. (7) yields 



TT , ^((25/6)-^) 
"^""^ ^ 2 + ln(A/c.) 



O 



1 



5ln (A/a>) 



(20) 



in the limit w — > 0. The bare value of g cancels out in 
the second term. Electron interactions provide therefore 
the leading universal logarithmic correction to scaling of 
conductivity at low frequencies in this case. Similar log- 
arithmic corrections arise in perturbative quantum chro- 
me dynamics, for example. 

To get an estimate of the size of the first-order cor- 
rection to conductivity in Eq. (19) we restore all the 
constants we previously set to unity. For the dielectric 
constant of e « 6 we find g ~ 2 The conductivity is 
determined by the bare value of disorder, not the inter- 
action, however. The disorder coupling can be written 
as A^/2 ^ $j/(I)o, where $^ w /if^ jg ^j^g fl^x of the ef- 
fective random magnetic field h — diA2 — d2Ai through 
the ripple of the average size ^, and $0 = hc/e is the fiux 
quantum. Effective magnetic field h can be estimated by 
assuming that the same random magnetic field is respon- 
sible for the observed suppression of weak localization in 
graphene 17|. Using /i « IT and ^ « 30r 

and fj « 3e^/h, to the leading order. 



17| we find 
Clearly 



A « 2, 

this is only a crude estimate and for A of order one the 
higher-order terms in Eq. (19) need to be included. It 
is encouraging, however, that the size of the lowest-order 
correction is significant and in the direction towards the 
experimental result. The latter result is a direct and 
a non-trivial consequence of the unscreened long-range 
Coulomb interaction in graphene. It may also be rele- 
vant that the experimental observation of normal local- 
ization properties in graphene seems to correlate with 
a lower minimal conductivity 18[. The present theory 
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would naturally account for this since both effects result 
from a low value of A. 

A limitation of our result should also be pointed out. 
We have computed the T = 0, ^ conductivity along 
the line of fixed points, whereas the measurements of 
the minimal dc conductivity in graphene typically corre- 
spond to the opposite, w = 0, T — > limit. While both 
conductivities at a fixed disorder A are expected to be 
universal, the two numbers could in principle be differ- 
ent [3]. If the results from the related critical bosonic 
theories with and without disorder [l^, Q are of 
any guide, the latter universal number may be expected 
to be only larger. Its computation along the line of the 
fixed points at which both the disorder and interactions 
are finite may be a non-trivial task though. 

Our result at A = is in stark contrast to the re- 
cent result in ref. [2l[, where Coulomb interaction is 
claimed to suppress the conductivity at low frequencies. 
We note that this conclusion is in contradiction to the 
well-established infrared irrelevancy of the Coulomb cou- 
pling [3| , 01 J which implies that the Coulomb interaction 
can only provide corrections to the Gaussian conductiv- 
ity, which ultimately vanish in the dc limit. Furthermore, 
the result of ref. 2]J is explicitly dependent on the cutoff 
in the Dirac theory, which is completely arbitrary. In con- 
trast, the cancelation of the two logarithmically divergent 
terms [22[ in Eq. (18) guarantees the renormalizability of 
our result. The low- frequency conductivity we computed 
is consequently perfectly cutoff-independent to the order 
of our calculation, as it has to be if the picture of Dirac 
quasiparticles is to remain physically meaningful in pres- 
ence of the interactions. 

Finally, let us adress the density of single-particle 



states: N{u;) 



OC UJ 



(2-z)lz 



just as at the unstable line 



at (7 = 0. Since g ^ however, Eq. (8) implies z = \ 
at the stable line. The density of states is thus exactly 
linear, in contrast to the 5 = line. 

To summarize, our main finding is that the lowest- 
order combined effect of electron interactions and rip- 
pling in graphene is to increase its minimal dc conduc- 
tivity in a non-universal, disorder-dependent fashion. A 
testable prediction of our theory would be a decrease of 
minimal conductivity in graphene with the suppression 



of wrinkling, which, incidentally, should also bring back 
the usual localization behavior at finite density. 
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